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ABSTRACT 

ARGO-YBJ is an air shower detector array with a fully covered layer of 
resistive plate chambers. It is operated with a high duty cycle and a large field 
of view. It continuously monitors the northern sky at energies above 0.3 TeV. In 
this paper, wc report a long-term monitoring of Mrk 421 over the period from 
2007 November to 2010 February. This source was observed by the satellite-borne 
experiments Rossi X-ray Timing Explorer and Swift in the X-ray band. Mrk 
421 was especially active in the first half of 2008. Many flares are observed in 
both X-ray and 7-ray bands simultaneously. The 7-ray flux observed by ARGO- 
YBJ has a clear correlation with the X-ray flux. No lag between the X-ray and 
7-ray photons longer than 1 day is found. The evolution of the spectral energy 
distribution is investigated by measuring spectral indices at four different flux 
levels. Hardening of the spectra is observed in both X-ray and 7-ray bands. 
The 7-ray flux increases quadratically with the simultaneously measured X-ray 
flux. All these observational results strongly favor the synchrotron self-Compton 
process as the underlying radiative mechanism. 

Subject headings: BL Lacertae objects: individual (Markarian 421) - gamma rays: 
observations 
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Introduction 



Mrk 421 {z = 0.031) is one of tlie brightest blazars known and is classified as a BL 



Lac object, a subclass of active galactic nuclei (AGNs). Mrk 421 was the first B 



source detected (by EGRET in 1991) at energies above 100 MeV (ILin et al 



Lac 



1992|), and 



was also the first extragala ctic obiect detecte d by a ground-based experiment (Whipple 



at energies around 1 TeV (jPunch et al. 



19921 ) (in the following we will refer to 7-rays as 



those around 1 TeV). Its emission, like that of the other blazars, is generally dominated 
by nonthermal radiation from a relativistic jet aligned along our line of sight. The 
spectral energy distribu tion (SEP) is doub le-humped at X-ray and 7-ray energies in 
a plot of uFi, versus u fiFossati et al.lll998l ). where u is the frequency and F^, the flux 
density. The hump at low energies is usually interpreted as being due to synchrotron 
radiation from relativistic electrons (and positrons) within the jet. The origin of the hump 
at high energies is under debate. Many models attribute the high-energy emission to 
the inverse Compton scattering of the synchrotron (synchrotron self-Compton, SSC) or 
external photons (external Compton , EC) by the same population of relativistic electrons 
f Ghisellini et al. 1998 : Dermer et al. 1992I ). therefore an X-ray/7-ray correlation would 
naturally be expected. Other models invoke hadronic processes including proton-initiated 
cascades and/or proton-synchrotron emission in a magnetic-field-dominated jet. Although 



( lAharonian 



2000 



Miicke et al. 


2003) 


fiGaidos et al 


1996) 



, they are generally challenged by the most rapid flares 



Mrk 421 is a very active b lazar with r naior outbursts about o nce every two years 



in both X-ray (Cui et al. 



20041 ) and 7-ray flTluczykont et al. 



2010) bands. A major 



outburst usually lasts several months and is accompanied by many rapid flares with 
timescales from tens of minutes to several days. Its high variability and broadband emission 
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require long-term, well-sampled, multiwavelength observations in order to understand 
the emission mechanisms of these outbursts. During the last decade, several coordinated 
multiwavelength campaigns focusing on Mrk 421 have been conduct ed both in response 



to strong out 



Fossati et al. 



jurst s and as part of dedicated observation campa i gns (iRebillot et al. 



20081; 



Acciari et al. 



2009 



Donnarumma et al 



2009 



Horan et al.l 



2006 



2009). Some 



important general features of the AGN flares have been obtained. Although X-rays and 
7-rays are found to be strongly correlated, neither type is evidently correlated with optical 



and radio em issions. The spectra l index becomes 



7-ray bands (IRebillot et al. 



2006 



larder at higher fluxes i n both X-ray and 



Krennrich et al 



2002 



Aielli et al. 



2010l ). An intensive 



multiwavelength monitoring campaign has recently been cond ucted with the Whipp 



20051). 



telescope and the Rossi X-Ray Timing Explorer (RXTE) (jBlazejowski et al. 
Similar features, including correlated variability at different energies, flaring and spectral 
evolution are also observed. All these phenomena can be interpreted in the framework of the 
SSC model. However, "orphan flares" , which have only 7-ray emission without low-energ y 



companions, and a lag of about two days between X-rays and 7-rays ( jBlazejowski et al. 
2OO5I ) are usually recognized as major challenges to the model. 



A long-term simultaneous X-ray/7-ray observation is better performed by means 



of a combination of satellite-borne X-ray experi ments and wide fle 



experiments, such as the Tibe t AS-7 experiment (lAmenomori et al. 



experiment ( 1 Aielli et al. 



d-of-view air shower 



20031) and ARGO-YBJ 



20061 ). which are operated day and night with a duty cycle higher 



than 85% and can observe any source with a zenith angle less than 50°. This is essential in 
order to investigate the temporal features of AGN emissions. The ARGO-YBJ experiment 
has continuously monitored the northern sky for outbursts from all AGNs, such as Mrk 421, 
since 2006 June. Meanwhile, these sources were also monitored by the satellite-borne X-ray 
detectors All-Sky Monitor (ASM) /RXTE and Burst Alert Telescope (BAT) /Swift. In this 
paper, we report on the long-term monitoring of Mrk 421 for 7-ray outbursts and on the 
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correlation between 7-rays and simultaneous X-rays over the period from 2007 November 
to 2010 February. The paper is organized as follows: the ARGO-YBJ experiment is briefly 
introduced in Section 2 and its long-term performance is shown in Section 3. A data 
analysis method is described in Section 4. Observation findings are presented in Section 5. 
Conclusions are given in Section 6. 



2. The ARGO-YBJ Experiment 

The ARGO-YBJ experiment, located in Tibet, China at an altitude of 4300 m a.s.l., is 
the result of a collaboration among Chinese and Italian institutions and is designed for very 
high energy 7-ray astronomy and cosmic ray observations. The detector consists of a single 
layer of resistive plate chambers (RPCs), which are organized with a modular configuration. 
The basic module is a cluster (5.7 m x 7.6 m) composed of 12 RPCs (2.850 m x 1.225 m 
each). The RPCs are equipped with pick-up strips (6.75 cm x 61.80 cm each), and the 
logical OR of the signal from eight neighboring strips constitutes a logical pixel (called a 
"pad" ) for triggering and timing purposes. One hundred thirty clusters are installed to form 
a carpet of about 5600 m^ with an active area of ~93%. This central carpet is surrounded 
by 23 additional clusters (a "guard ring") to improve the reconstruction of the shower core 
location. The total area of the array is 110 m x 10 m. More details about the detector 



and RPC performance can be found in, for example, lAielli et al. 



fl2006h 



The RPC carpet is connected to two independent data acquisition systems 
corresponding to two different operation modes, referred to as the shower and the scaler 



( lAielli et al. 



20081 ) modes. Data used in this paper refer to the shower mode, in which 
the ARGO-YBJ detector is triggered when at least 20 pads in the entire carpet detector 
are registered within 420 ns. The high granularity of the apparatus permits a detailed 
spatial— temporal reconstruction of the shower profile and therefore the incident direction 
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of the primary particle. The arrival time of the particles is measured by time to digital 
converters (TDCs) with a resolution of a pproximately 1 .8 ns. In order to calibrate the 



18,360 TDC channels, an off-line method (iHe et al. 



20071 ) has been developed using cosmic 



ray showers. The c alibration precision is 0.4 ns, and the procedure is applied every month 



( lAielli et al. 



2009ah . 



The central 130 clusters began taking data in 2006 June, and the "guard ring" was 
merged into the DAQ stream in 2007 November. The trigger rate is ~3.6 kHz with a dead 
time of 4%, and the average duty cycle is higher than 85%. 



Detector Performance 



For long-term monitoring campaigns, the stable operation of the equipment is very 
important. In order to continuously monitor the performance of the RPCs, including 
detection efficiency and time resolution, a cosmic ray muon telescope is set up near the 
detector array. The RFC efficiency fluctuates by about 0.3% and the time resolution by 
about 0.4 ns in a day, and these values become 1.5% and 1 ns in a year, respectively. 



Detailed in: 



Aielli et al 



brmatio n about the performance monitored using this telescope can be found in 



fl2009ch . 



To estimate the angular resolution and effective full Monte Carlo simulation o 



the RFC detector array is developed. In the code, the C ORSIKA packag e flHeck et al 



i s used to describe the a ir shower development. G4argo (iGuo et al. 



19981 ) 



( lAgostinelli et al. 



20101 ) ■ a GEANT4-based 



20031 ) package, is used to simulate the response of the RFC array. For 
events with a number of fired pads (Npad) greater than 100, the Foint Spread Function 
(FSF) has a single Gaussian functional form. For events at lower Np^^, the best fit to the 



FSF becomes a combination of two Gaussian distributions, the wider of which contains 
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20% of the events. To simplify the description of the PSF, a parameter tpjo is defined as 
is the opening angle containing 71.5% of the events. When the PSF is a single Gaussian, 
ipYO maximizes the signal-to-background ratio for a point source. For Npad > 1000, ipjo is 
0.47°, while at Npad ~ 20 ipio becomes 2.8°. The effective area of the detector for 7-induced 
showers depends on the 7-ray energy and incident zenith angle, e.g., it is about 100 m^ at 



100 GeV and >10,000 m^ above 1 TeV for a zenith angle of 20° flAielli et al. 



2009bf ). 



The angular resolution, pointing accuracy and stability of the ARGO-YBJ detector 
array have been th oroughly tested by measuring the shadow of the Moon in cosmic rays 



( lluppa et al 



2OO9I ). The shadow is detected with a significance of 10 a per month using 



the ARGO-YBJ data. The position of the shadow allows the investigation of any pointing 
bias. The east-west displacement is in good agreement with the expectation, while a 0.2° 
pointing error toward the north is observed and is under investigation. 



4. Data Analysis 

For the analysis presented in this paper, only events with a zenith angle less than 45° 
are used, and the data set is divided into six groups according to Npad- The event selections 
are listed in Table 1, where R is the distance between shower core position and the carpet 
ce nter, and TS is the time spread of the shower front in the conical fit defined in Eqation(l) 



of 



Aielli et al. 



( l2009al ). With these selections, the angular resolution is improved, e.g., for 
events with Npad > 60 and Npad > 100, the opening angle ipYo decreases from 1.68° and 
1.27° to 1.36° and 0.99°. As a consequence, the significance of the Crab Nebula is increased 
by about 10% and 25%, respectively. 

In order to obtain a sky map using events in each Npad group, an area centered at the 
source location in celestial coordinates (right ascension and declination) is divided into a 
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grid of 0.1° X 0.1° bins and filled with detected events according to their reconstructed 
origin. The number of events in each grid bin is denoted as n^, where the subscript i denotes 
the bin number. In order to e xtract an excess of 7-rays from the source, the direct integral 



method ( Fleysher et al. 



20041 ) is applied to estimate the number of cosmic ray background 
events in the bin, denoted as bi. An essential assumption in this estimation is that the 
background must be uniform around the source. However, an ani sotropy of the cosmic ra y 



flux is measured o ver spatial scales such as 10° x 10° and larger (lAmenomori et al 



Zhang et al. 



20061 : 



20091 ). This anisotropy as measured by the rii/bi ratio is stable; therefore, it 
is possible to correct it with a long-term measurement for each grid bin. An average of 
the ratio over the bins in a window 11° x 11° centered on the source bin is applied for 
smoothing. In this procedure, in order to avoid any contamination of the excess in the 
source bin and possible spread out due to the finite angular resolution, the contribution 
from a 5° X 5° window around the source bin is excluded. Finally the correction factor, 
denoted as is calculated as follows: 



^ m 



Tin 



m ^ — ' hi 

j=l J 



where the subscript j is the index of the m = 12100 — 2500 = 9600 selected grid bins. The 
corrected number of background events in each bin is b* = /3ibi. The typical value of /3 
around Mrk 421 is approximately 0.9995. The value of /3 for each bin is calculated using 
about two years of data and is stored in a database for routine analysis. 

Taking into account the PSF of the ARGO-YBJ detector, the events in a circular area 
centered on the bin with an angular radius of ipjQ are summed together. Namely, 



i=l 



1=1 



where k is the number of bins in the circular area, Non is the total number o: 



Nf) is the number of background events. The Li— Ma formula (ILi &: Ma 
estimate the significance. 



(2) 



events, and 



19831 ) is used to 
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5. Results 

The data used in this paper were collected by the ARGO-YBJ experiment in the 
period from 2007 November to 2010 February. The total lifetime is 676.0 days. The 
numbers of events in different groups after the selections are listed in Table 1. A clear signal 
from Mrk 421 with significance greater than llcr is observed using events with Npad > 60 
(see Figured]). A signal at such a level of significance allows us to study flux variations, 
correlations with the X-ray flux, and the evolution of the SED. 



5.1. Temporal Analysis 

In order to study the correlation between 7-rays and X-rays, the daily averaged light 
curves of both the hard X-rays (15—50 keV) measured by BAT/Swz/jl and the soft X-rays 
(2—12 keV) measured by ASM/i?XT£cl are used. The observations by RXTE and Swift 
have a rather long exposure by orbiting the Earth every 1.5 hr. Since the fluctuation of the 
X-ray flux is abnormally large in some days, in order to control the quality of the data, days 
that have a very large error on the mean daily event rate are removed from the data set. For 
ASM/ RXTE, the distribution of the error indicates that a selection of the errors smaller 
than 1 count s^^ will cut everything beyond four standard deviations in the distribution. 
A similar cut applies to the BAT/ Swift data, in which a selection of the errors smaller 
than 0.0035 counts cm~^ s~^ cuts everything beyond four standard deviations in the error 
distribution. Approximately, 6.4% and 5.6% of events are removed from the RXTE and 



^ Transient monitor results provided by the BAT/ Swift team: 
http : / /heasarc . gsf c . nasa . gov/docs/ swift/ results/transients/weak/Mrk42T7J 

^Quick-look results provided by the ASM/ RXTE team: 
http : //xte -mit . edu/ASM_lc . htmlj 
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Swift data sets, respectively. Whether it is day or night, ARGO-YBJ observes Mrk-421 
while the AGN is in its field of view. A typical transit lasts usually 6 hr. An observational 
time less than 5 hr day~^ indicates some malfunctioning of the detector in that day, which 
is thus removed from the data set. In total, 9.7% of data are removed in this way. Finally, 
737, 728, and 712 days are selected from the ASM, BAT, and ARGO-YBJ reconstructed 
data sets, respectively. 



5.1.1. Light Curves 



In 552 days all three experiments observed Mrk 421 simultaneously. In Figure [2], the 
accumulation of event rates from the Mrk 421 direction is shown. The Swift event rate has 
been normalized using the RXTE scale and the ARGO-YBJ curve is obtained using events 
with Npad > 100, thus the median energy of the observed photons is 1.8 TeV, assuming a 
spectral index —2.4. The fast increase in the three curves indicates that the source had a 
long-term outburst at the beginning of 2008. The following quiet state lasted for about 
200 days. Afterwar d Mrk 421 became increasingly more active. In fact, there were fiares 



in 2009 November (j Isobe et al. 



2010). The duty cycle of ARGO-YBJ was low due to 



detector maintenance, there: 



2010 February (I Isobe et al 



ore i t is not obv ious in Figure [21 There was a large fiare in 



2010 



Ong 



2010|) 



Out of the long-term variation that is clearly revealed in the cumulative light curve 
shown in Figure [21 Mrk 421 undergoes a large outburst during the period from 2008 
February to June, indicated by the steepest part of the curves. In fact, it is a combination 
of several large fiares. A better view of these is shown in Figure [3l where a smoothing 
analysis is applied for both 7-ray and X-ray curves, and each point is the event rate 
averaged over five days. Four large fiares are observed by all three detectors, and the peak 
times are in good agreement with one another. The fourth fiare has been reported by the 
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ARGO-YBJ experiment in 



Aielli et al. 



teOlOl ). It gives an important observation when the 



Cherenkov telescopes are hampered by the Moon. It can be concluded that there exists a 
good long-term correlation between 7-rays and X-rays. 



5.1.2. X-ray /TeV Correlation 



The discrete correlation function (DCF) (lEdelson fc Krolidll988l ) is used to quantify 
the degree of correlation and the phase differences (lags) in the variations between 7-rays 
and X-rays. The daily fluxes before smoothing are used for this analysis. The DCF (in 
1 day bins) derived from RXTE and ARGO-YBJ data (with Npad > 100) is shown in 
the left panel of Figure HI where a positive value means that 7-rays lag X-rays. The peak 
of the distribution is around zero and the correlation coefficient at zero is ~ 0.77. The 
result derived from Swift and ARGO-YBJ data is shown in the right panel of Figure H] 
and the correlation coefficient at zero is ^0.78. To e stima te the lag and its uncertainty, 
a data-based simulation suggested by 



Peterson et al. 



(119981 ) is applied and the correlation 



coefficient between —10 and 10 days is fitted with a Gaussian function. The median 
value and corresponding 68% confidence level errors are — 0.14lQ g5 and — 0.941}'q7 days 
for the correlations of ARGO-YBJ/i^XTE and ARGO-YBJ / Swift data, respectively. No 
significant lag longer than one day is found. 



5.2. Spectral Energy Distribution 

To study the SED at different flux levels, the data simultaneously observed in 7-ray 
and X-ray bands are divided into four groups according to the observational time periods 
in which the ASM/ RXTE counting rate is — 2, 2 — 3, 3 — 5 or > 5 cm~^ s~^. For each 
group, a flux-averaged SED is constructed both at 7-ray and X-ray energies. 



- 14 - 



5.2.1. X-ray Spectra 



ASM/ RXTE monitors the X-ray emission fro m Mrk 421 at three energy bands, i.e. 



1.5 — 3, 3 — 5 and 5 — 12 keV (iLevine et al. 
column density 1.38 x 10^° cm~^ (jDickey &: Lockman 



19961). In the flux estimation, the hydrogen 



19901 ) and a power law spectrum are 



assumed. The best-fit spectral indices for the four flux levels are —2.43 ±0.04, —2.15 ±0.03, 
—2.05 ± 0.03, and —2.02 ± 0.08, respectively, in which only statistical errors are taken 
into account. The spectral indices versus the correspondi ng fluences at 1 keV are shown 



Rebillot et al. 



toom . in which 



in Figure |5l This result is consistent with the analysis of 
a spectral hardening toward high fluxes is also reported based on a shorter timescale 
observation. This indicates that this correlation is independent of the timescale. 



5.2.2. ^-ray Spectra 



To estimate the spectrum of 7-rays with a distribution of the number of events in 



excess as a function of iVpn^, we fo ^ 



elsewhere (lAmenomori et al. 



2009 



low a wide^ 



Aielli et al. 



used method that is described in detail 



2010). In this procedure, we assume for the 



spectrum of Mrk 421 a power law with a cutoff factor e which takes into account the 
absorption of 7-rays in the extra g alacti c background light. We adopt the optical depth t(E) 



estimated by 



Franceschini et al. 



(120081 ). The ARGO-YBJ detector response is also taken 
into account. The simulated events are sampled in the energy range from 10 GeV to 100 
TeV. 

To test this method, the same analysis is performed with the data in the direction 
of the Crab Nebula, the standard candle in the 7-ray sky. The resulting spectrum is 
(4.2 ±0.4^tot) X lO-^HE/TeV)"^- ^^^"-"^"'"* photons TeV"^ cm"^ g-i^ ^j^jj.^ jg agreement 



with our previous measurement (lAielli et al. 



2OIOI ) and observations by other detectors. 
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(lAmenomori et al. 



20091). 



such as H.E.S.S. flAharonian et all 120061). MAGIC flAlbert et all 120081 ). and Tibet AS-7 



Applying this procedure to Mrk 421, we obtain the spectra for the four event groups 
with different flux levels. The spectral indices in the energy range from 300 GeV to 10 TeV 
are -2.48 ±0.22, -2.53 ±0.21, -2.15 ±0.18, and -1.87±0.21, respectively. Only statistical 
error is quoted. The corresponding flux above 1 TeV ranges from 0.8 to 6 times that of the 
Crab Nebula unit, i.e., 2.67 x 10^^^ photons cm~^ s~^. The spectra seem to become harder 



with increasing flux, as indic ated in Figure [6 



by the W hipple experiment (IKrennrich et al. 



elsewhere ( lAielli et al. 



in a greement with the function obtained 



2OO2I ) ■ A similar result has been reported 



I2OIOI ) using the three-day flare data in 2008 June. The quoted errors 
in Figure [6] are statistical. The s ystematic error is estimated to be <30% in the flux level 



determination (lAielli et al. 



2010) 



5.2.3. Correlation Between ^-ray and X-ray Fluxes 



Using the spectra described above, we investigate the correlation between 7-ray and 
X-ray fluxes. Figure [7] shows the integral 7-ray flux above 1 TeV as a function of the 
integral X-ray flux from 2 keV to 12 keV; a positive correlation is observed. A quadratic 
flt (with the function y = ax^ ± 6 ) to the data points yields x^/dof= 1.9/2, while a linear 
flt yields x^/dof= 7.7/2, where dof refers to degrees of freedom. The observation favors a 
quadratic correla tion between 7 -ray a nd X-ray fluxes. A similar quadratic correlation has 



been report ed bv iFossati et al. 



obtained by 



Amenomori et al 



2OO8I) . In contrast, a n observation with 



(120031). According to 



Katarzyhski et al. 



inear correlation is 



(120051), changes of 



the magnetic fleld, electron density, and adiabatic cooling may be associated with different 
correlations between 7-ray and X-ray fluxes. 
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5.3. Modeling of the X-ray and 7-ray Emissions 



A fit to the four flux-averag e d SED s with a 



prop osed by 



homogeneous one-zone SSC model 



Mastichiadis &: Kirk 



1997 



Yang et al. 



Mastichiadis fc KirkI (119951 ) (see also 
20081 ) is performed. In this model the parameters include the Doppler factor 5 = 
l/[r(l — f3cos9)]; the spherical blob radius R; magnetic fleld strength B; electron spectral 
index s; electron maximum Lorentz factor 7max; and electron injection compactness 
le = ^meCa-^R^ did ~ ^)Qe, where T and cf3 are the Lorentz factor and the speed of 
the blob, respectively, 7 is the electron Lorentz factor, (Tt is the Thomson cross section, 6 is 
the angle between its direction of motion and the line of sight of the observer, and Qe, the 
electron spectrum at injection, is assumed to follow a power law Qe = ?e7~'' exp(— 7/7inax)- 
The best flts are shown in Figure [8] for different flux levels, with the corresponding 
parameters given in Table 2. In our flts, t he magnetic fleld str e ngth is estimated by 



Mastichiadis fc KirkI (119971 )). where i/g.is is the 



B = 5 X 10^^5z/s,i8'^(r27 (s^^ Equation (6) of 
synchrotron peak frequency in units of 10^^ Hz and z^c,27 is the IC peak frequency in units of 
10^'' Hz. For the lowest flux level (see Figure (Ht^l)), the magnetic fleld strength is estimated 
to be ~ 0.08 G. Compared to the lowest flux level. X-ray peak frequencies in other flux 
levels increase by a factor of ~ 2, but the IC peak frequencies have few changes. Therefore, 
the magnetic fleld strengths in other flux levels are larger than that in the lowest flux level 
by a factor of ~ 2 when the Doppler factor 6 is roughly flxed. 



6. Discussion And Summary 

Mrk 421 is a very active blazar with frequent outbursts, which are composed of many 
flares and can last as long as a few months. This makes this blazar an excellent candidate 
for studying the jet physics in AGNs. A strong correlation between its 7-ray and X-ray 
emissions has been conflrmed by many observations in the past decade (for a review see 
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Wagner I (120081 )). Most of the previous 7-ray observations, however, are carried out by 
Cherenkov telescopes with hmited exposure and usually focus on short timescales. In 
contrast, the high duty cycle of the ARGO-YBJ experiment makes possible a long-term 
and continuous observation of this variable source, allowing simultaneous monitoring of 
7-rays and X-rays for about two years. This increases the set of long-term simultaneous 
multiwavelength observations of Mrk 421, which are essential for studying the correlation 
between energy bands where different emission mechanisms are at work. The observation 
time, from 2007 December to 2010 February, covers both active and quiet phases. The 
7-ray flux shows a good long-term correlation with the X-ray flux (see Figure and all the 
large X-ray flares have their 7-ray counterparts during the outburst time (see Figure |3]), 
i ndicating that 7-rays and X-rays may have a common origin as assumed in the SSC model 



(iMastichiadis fc Kirk 



I995I). 



In the SSC model, the 7-ray photons are produced via inverse Compton scattering 
off the synchrotron photons by the same electrons, and simultaneous variability or short 
lags are expected between 7-ray and X-ray fluxes. Short lags can be caused by differences 
in acceleration and cooling timescales or by reverse shocks, and sub-hour la gs have been 



definitely measured between different X-ray (IRavasio et al. 



20071 ) interband energies and between X-ray and 7-ray bands ( iFossati et al. 



20041) and ^-rav 



(lAlbert et al. 



2008 



). On the 



other hand, the characteristic timescale of the SSC process would be to o short to account 



Blazejowski et al 



( 120051 ) with a marginal 



for a lag of two days such as that reported by 
significance. In this paper, the two-year data are used to search for possible lags between 
variations in X-rays and 7-rays. No lag longer than one day is observed (see Figure S]). 

A sudden variation of fiux can be caused by different reasons, e.g., a change in the 
number of emitting electrons and/or the maximum momentum of emitting electron and/or 
the magnetic field strength, with different evolutions of the SED at X-rays and 7-rays in the 
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SSC model. To investigate the evolution at different ffux levels, both the 7-ray and X-ray 
data are divided into four groups according to the X-ray flux. A hardening in the spectra 
toward high fluxes is observed (see Figure O and Figure (HI). The results base d on two years of 



data are consistent with the results obtained by the Whipple experiment (IKrennrich et al. 



20021 ). A close tie between the variation of the flux and of the sp ectral index indicate s 



peak energy increases with the flux, which has also been found in lAleksic et al. 



mm . 



This supports the pre diction of the SSC model 



of emitting electrons (IMastichiadis &: Kirk 



br changes in the maximum momentum 



19971 ). Moreover, we study the correlation 



function between 7-rays and X-rays, finding that the 7-ray flux shows a quadratic increase 
with the X-ray flux (see Figure [7]). In the homogeneous SSC model, the synchrotron flux 
is proportional to the electron density, and the IC 7-ray flux is proportional to both the 
electron density and the synchrotron flux; therefore, the 7-ray flux is a quadratic function 
of the synchrotron flux. Never before was there an indication distinguishing quadratic from 
linear correl ations between 7-ray and X-ray fluxes according to observations of flares, as 



reviewed by 



Wagner 



( I2OO8I ). We also construct a homogeneous one- zone SSC model to 
simultaneously fit the 7-ray and X-ray emissions in four different flux levels (see Figure [H]) 
by changing the electron parameters Ig and/or 'jmax and/or magnetic field strength. We 
find that the flux variation seems to be caused by the variation of the maximum energy and 
density of the electron injection spectrum. 

In conclusion, we have presented a long-term continuous monitoring of Mrk 421 and 
a correlation between 7-rays observed by the ARGO-YBJ experiment and satellite-borne 
X-ray data. The temporal and spectral analysis strongly support the predictions of the SSC 
model. 
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Table 1. Event Selections and the Number of Events 





/? fml 

± Kj I J-J-J- I 


-L kJ I llO 1 


iMninhpT* of r AT'P'n'l"'^; 


[20, 60] 


No cut 


<50 


8.71x10^° 


[60, 100] 


No cut 


<30 


1.83x10^° 


[100, 200] 


R <70 


<20 


6.13x10^ 


[200, 500] 


R <70 


<20 


3.42 xlO'^ 


[500, 1000] 


R <60 


<20 


1.05 xlO*^ 


[> 1000] 


R <30 


<20 


3.71x10^ 



Table 2. Best-Fit Parameters in the SSC Model 



Flux Level ^max U B (G) R (cm) 5 a 



1 


7 X 10^ 


6 X 10-6 


0.08 


5 X 10^6 


16 


1.7 


2 


7 X 10^ 


1 X 10-5 


0.15 


5 X 10^6 


15 


1.7 


3 


1 X 10^ 


1 X 10-5 


0.15 


5 X 10^6 


15 


1.7 


4 


2 X 10^ 


1.4 X 10-5 


0.15 


5 X 10^6 


15 


1.7 
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Fig. 1. — Distribution of statistical significance around Mrk 421. 
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Fig. 2. — Cumulative light curves from the Mrk 421 direction. The red curve is the 7-ray 
result observed by ARGO-YBJ, and the shaded red region indicates the corresponding la 
statistical error; the black curve represents soft X-rays (2—12 keV) observed by ASM/ RXTE. 
Hard X-rays (15—50 keV) observed by BAT/ Swift are given by the blue curve, and the scale 
has been normalized to the ASM/ RXTE one. 
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Fig. 3. — Daily light curves from Mrk 421 direction in different energy bands from 2008 

February 1 to July 18. Each bin contains the event rate averaged over the five-day interval 
centered on that bin. The panels from top to bottom refer to 2—12 keV [ASM /RXTE), 
15-50 keV (BAT /Swift), and 7-ray (ARGO-YBJ), respectively. 
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Lag [Day] Lag [Day] 

Fig. 4. — Discrete correlation function between X-ray and 7-ray light curves from 2007 
November to 2010 February. Left: 2-12 keV {ASM/ RXTE) vs. 7-ray (ARGO-YBJ); a 
Gaussian function is used to fit from —10 to 10 days. Right: 15—50 keV (BAT /Swift) vs. 
7-ray (ARGO-YBJ). Positive value means that 7-rays lag X-rays. 
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Fig. 5. — Correlation between the X-ray flux at 10 keV and the corresponding photon index 
at 2-12 keV. 
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Fig. 6. — Spectral index vs. 7-rav flux above 1 TeV . The solid line is the function obtained 
by the Whipple experiment (IKrennrich et al.ll2002l ). 



-29- 




Fig. 7. — 7-ray flux above 1 TeV vs. X-ray flux at 2—12 keV. The solid line is a quadratic 
fit using function y = ax^ + 6, which yields x^/dof= 1.9/2. The dotted line is a linear fit, 
which yields x^/dof=: 7.7/2, where dof refers to degrees of freedom. 
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Fig. 8. — Spectral energy distribution of Mrk 421. (1) to (4) are derived from four flux level 
data groups from low to high according to the ASM/RXTE counting rate (see the text for 
details). The solid line shows the best fit to the data with a homogeneous one-zone SSC 
model, and the best- fit parameters are hsted in Table 2. 



